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In  the  area  of  aircraft  structural  dynamic  design,  the  existence  of  a  reliable  structural 
optimization  method  can  be  a  valuable  design  tool,  enabling  a  designer  to  perform  rapid 
identification  and  evaluation  of  the  most  promising  structural  elements,  which,  upon 
modification  of  their  modal  properties,  could  minimize  vibration  response.  In  this  pro¬ 
gram,  a  feasibility  study  was  performed  to  evaluate  two  structural  optimization  methods, 
the  Vincent  circle  and  the  forced  response  strain  energy  method  for  reducing  vibration 
response,  primarily  via  structural  stiffness  changes,  using  NASTRAN  beam-element  repre¬ 
sentation  of  the  AH-1G  with  uniform  damping.  The  forcing  function  consisted  of  a 
single  main  rotor  2  rev  steady-state  vertical  excitation.  The  pilot's  seat  was  selected  as 
the  location  for  dynamic  response  reduction. 

The  relative  effectiveness  of  the  methods  to  minimize  the  dynamic  response  at  the 
selected  location,  subject  to  the  constraints  and  considerations  stated  above,  was  based 
essentially  on  two  factors:  (a)  selection  of  practical  stiffness  changes  that  an  element 
can  undergo,  and  (b)  evaluation  of  the  effectiveness  of  the  groups  of  elements  predicted 
by  each  method  by  calculating  the  response  of  the  structure  subject  to  practical  stiffness 
changes. 

The  results  obtained  from  the  two  methods,  i.e.,  stiffness  changes  required  and  group 
of  elements  identified  with  high  potential  for  reducing  vibration,  at  best,  do  not  agree. 
These  results  indicate  that  the  forced  response  strain  energy  method  is  more  suitable  in 
attaining  response  reduction  according  to  the  investigative  methodology  used  in  this 
program. 

The  Vincent  circle  method  was  further  investigated  in  conjunction  with  mass  and  damp¬ 
ing  changes.  Also,  the  possibility  that  the  Vincent  circle  method  would  indicate  the 
optimum  location  of  a  fixed  mass  and  damping  dynamic  absorber  was  investigated, 
yielding  negative  results. 

The  project  engineer  for  this  effort  was  Mr.  N.  J.  Calapodas  of  the  Aeronautical  Tech¬ 
nology  Division,  Structures  Technical  Area.  , 
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model  of  the  AH-1G  helicopter.  The  forced  response  strain  energy 
method  was  then  applied  to  a  large  complex  built-up  NASTRAN  AH-1G 
model.  This  application  provided  useful  comparative  data 
identifying  the  structural  elements  considered  to  be  the  primary 
contributors  to  the  response.  Realistic  structural  stiffness 
changes  in  these  elements  were  assessed  to  determine  their 
effect  on  vibration  reduction.  The  Vincent  circle  method  was 
further  evaluated  for  mass  tuning,  damping,  and  dynamic  absorber 
parameters  using  the  elastic-line  model. 
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PREFACE 


This  investigation  of  vibration  reduction  through  structural 
optimization  was  performed  under  Contract  DAAK51-78-C-0011  for 
the  Applied  Technology  Laboratory,  U.  S.  Army  Research  and 
Technology  Laboratories  (AVRADCOM),  Fort  Eustis,  Virginia. 

The  program  was  implemented  under  the  technical  direction  of 
Mr.  Nicholas  J.  Calapodas  of  the  Applied  Technology  Laboratory. 

A  literature  survey  of  structural  optimization  techniques  for 
vibration  reduction  was  performed  and  two  techniques  -  the 
Vincent  circle  method  andtfae  forced,  response  strain  energy 
method  -  were  evaluated/tlsTng  a  NASTRAN  analysis  of  the  Model 
AH-1G  helicopter.  Mr.  (Horace  W.  Hanson  was  the  Bell  Heli¬ 
copter  Textron  Project  Engineer.  _ _ 
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INTRODUCTION 


In  the  development  of  a  helicopter  airframe  structure,  part  of 
the  design  cycle  must  be  concerned  with  keeping  vibration  to  a 
minimum.  The  analysis  and  control  of  helicopter  airframe 
dynamic  characteristics  are  quite  complex.  The  reasons  are, 
in  part: 

The  large  number  of  main  rotor  harmonic  excitations 
simultaneously  present,  as  well  as  excitations  from 
secondary  sources  such  as  the  tail  rotor  and  shafting. 

The  large  variation  in  main  rotor  oscillatory  hub  forces 
and  moments  and  phase  relationships  as  a  function  of  ship 
gross  weight  and  center-of-gravity ,  airspeed,  temperature, 
altitude,  etc. 

The  large  number  of  locations  at  which  vibrations  must  be 
controlled  within  the  aircraft  for  reasons  of  crew  and 
passenger  comfort,  system  reliability,  and  service  life 
of  the  airframe  structure  and  components. 

Design  changes  for  vibration  control  that  must  be  accom¬ 
plished  within  the  framework  of  the  overall  design  cri¬ 
teria,  including  many  considerations  such  as  size, 
strength,  weight,  and  aerodynamic  drag. 

Vibration  reduction  is  primarily  accomplished  through  dynamic 
devices  such  as  vibration  isolation  systems  and  dynamic  ab¬ 
sorbers,  and/or  through  structural  modifications  that  improve 
the  vibratory  elastic  deformation  characteristics  of  the  air¬ 
frame.  Rapid  and  efficient  techniques  are  needed  for  evalua¬ 
ting  and  optimizing  these  structural  changes  to  minimize  the 
vibration  response.  Two  methods  of  reducing  vibration  through 
structural  optimization  were  evaluated  in  this  study  -  the 
Vincent  circle  method  and  the  forced  response  strain  energy 
method. 
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LITERATURE  SURVEY 


A  brief  search  was  conducted  to  locate  and  review  existing 
literature  applicable  to  vibration  reduction  by  structural 
optimization  techniques  with  emphasis  on  the  Vincent  circle 
and  the  forced  response  strain  energy  methods.  Three  data 
bases  were  '  sed  in  the  survey: 

-  The  National  Technical  Information  Service  (NTIS) 

The  Defense  Technical  Information  Center  (DTIC) 

The  Engineering  Index  (COMPENDEX) 

The  information  that  is  cataloged  in  the  National  Technical 
Information  Service  is  derived  from  publications  of  the  fed¬ 
eral,  state,  and  local  government  agencies;  private  industry; 
and  universities,  with  emphasis  placed  on  commercial  applica¬ 
tions.  The  information  available  from  DTIC  is  obtained  from 
many  of  the  same  sources  as  NTIS;  however,  the  emphasis  is 
placed  on  military  and  defense  usage.  The  Engineering  Index 
contains  publications  of  the  engineering  societies,  such  as 
the  proceedings  of  conferences,  journals,  and  magazines. 

Five  documents  applicable  to  the  Vincent  circle  method  were 
found  (References  1  through  5).  Only  one  report  addressing 


1 Vincent ,  A.  H.,  A  NOTE  ON  THE  PROPERTIES  OF  THE  VARIATION 
OF  STRUCTURAL  RESPONSE  WITH  RESPECT  TO  A  SINGLE  STRUCTURAL 
PARAMETER  WHEN  PLOTTED  IN  THE  COMPLEX  PLANE,  Westland  Heli¬ 
copters  Limited  Report  GEN/DYN/RES/010R,  Yeovil,  Somerset, 
England,  September  1973. 

2 Done,  G.  T.  S.,  and  Hughes.  A.  D. ,  THE  RESPONSE  OF  A  VIBRA¬ 
TING  STRUCTURE  AS  A  FUNCTION  OF  STRUCTURAL  PARAMETERS, 
Journal  of  Sound  and  Vibration,  Vol .  38,  No.  2,  1975, 
pp.  255-266. 

3 Done ,  G.  T.  S.,  and  Hughes,  A.  D. ,  REDUCING  VIBRATION  BY 
STRUCTURAL  MODIFICATION,  Vertica,  Vol.  1,  No.  1,  1976, 
pp.  31-38,  (paper  presented  at  the  First  European  Rotor- 
craft  and  Powered  Lift  Aircraft  Forum  at  Southampton, 
September  22-24,  1975). 

4 Done,  G.  T.  S.,  Hughes,  A.  D.,  and  Webby,  J. ,  THE  RESPONSE 
OF  A  VIBRATING  STRUCTURE  AS  A  FUNCTION  OF  STRUCTURAL  PARA¬ 
METERS  -  APPLICATION  AND  EXPERIMENT,  Journal  of  Sound  and 
Vibration,  Vol.  49,  No.  2,  1976,  pp.  149-159.  . 

5Balmford,  D.  E.  H.,  THE  CONTROL  OF  VIBRATION  IN  HELICOPTERS, 
Aeronautical  Journal,  Vol.  81,  No.  794,  February  1977, 
pp.  63-67. 
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the  forced  response  strain  energy  method  was  located  (Refer¬ 
ence  6 ) . 

The  search  did  not  disclose  any  other  literature  that  dealt 
specifically  with  vibration  reduction  through  structural 
optimization. 

6Sciarra,  J.  J.,  USE  OF  THE  FINITE  ELEMENT  DAMPED  FORCED 
RESPONSE  STRAIN  ENERGY  DISTRIBUTION  FOR  VIBRATION  REDUCTION, 
Boeing-Vertol  Company  Report  D210-10819-1 ,  U.S.  Army  Re¬ 
search  Office  -  Durham,  Durham,  North  Carolina,  July,  1974. 
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STRUCTURAL  OPTIMIZATION  METHODS  UNDER  EVALUATION 


VINCENT  CIRCLE  METHOD 

If  a  damped  linear  structure  is  excited  by  a  constant  sinu¬ 
soidal  force  at  a  single  location,  the  response  of  some  other 
point  in  the  structure  can  be  shown  to  trace  out  a  circular 
locus  in  the  cqmplex  plane  when  any  single  structural  element 
stiffness  or  mass  parameter  is  continuously  varied  from  minus 
infinity  to  plus  infinity,  as  illustrated  in  Figure  1.  This 
dynamic  property  of  linear  structures  is  known  as  the  Vincent 
circle  phenomenon  (Reference  1).  The  phenomenon  requires  that 
some  value  of  damping  be  present  in  order  to  develop  the 
imaginary  component  of  the  response.  The  minimum  response 
attainable,  due  to  a  particular  stiffness  or  mass  parameter, 
is  indicated  by  the  radius  of  the  response  circle  and  location 
of  its  center.  Use  of  these  circular  properties  will  here¬ 
after  be  referred  to  as  the  Vincent  circle  method. 

The  relevant  mathematics  of  this  circular  response  phenomenon, 
based  on  the  theory  originally  established  by  Vincent,  are 
given  by  Done  and  Hughes  (Reference  2)  and  are  summarized  as 
follows : 

The  structure  shown  in  Figure  1  has  many  degrees-of-freedom, 
and  the  response  at  point  q  due  to  a  single  sinusoidal  force 
at  point  p  is  being  examined.  The  amplitude  of  the  displace¬ 
ment  response  vector,  x,  for  the  structure  is  given  by 

x  =  GF  (1) 

where 

"  2  _1 

G  =  [K  -  Muj  +  ituC]  (2) 

Consider  two  points,  r  and  s,  in  the  structure  having  mutually 
compatible  degrees  of  freedom.  A  linear  spring  of  stiffness, 
k,  is  inserted  between  these  two  points  so  as  to  exert  zero 
force  when  the  original  structure  is  in  equilibrium.  Con¬ 
sidering  the  structure  as  a  free  body,  the  forces  exerted  on 
it  at  points  r  and  s  by  the  spring  are  F  and  F  ,  respectively, 
which  have  the  relationship  r  s 

Fr  =  k<*s  "  V  =  -Fs  O) 
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Apply  constant  sinusoidal  force  of  amplitude  F  at 
point  p,  vary  spring  stiffness  k  between  points 
r  and  s,  and  measure  response  x  at  point  q. 


Imaginary 


Figure  1.  Vincent  circle  phenomenon. 


The  nonzero  elements  in  the  forcing  vector  F  are  F  ,  F  , 
F  .  The  response  equations  of  interest  are 

x  =  G  F+G  F+G  F 
q  qp  p  qr  r  qs  s 


X  =  G  F+G  F+G  F 
r  rp  p  rr  r  rs  s 


X—  G  F+G  F+G  F 
s  sp  p  sr  r  ss  s 


where 


G^  is  the  complex  receptance  providing  the  dis¬ 
placement  at  point  q  due  to  a  force  at  point  p, 
etc . 


Expressing  Fr  and  Fg  in  terms  of  xr  and  xg  from  Equation  (3), 
and  subsequent  elimination  of  xr  and  xg  from  Equations  (4), 

( 5 ) ,  and  ( 6 ) ,  gives 


J  =  G  +  H  (GsP  -  Grp)  (Gqr  -  Gqs) 

FP  ®  1  +  k<Grr  +  Gss  -  Grs  *  Gsr> 

This  may  now  be  rewritten  in  more  general  terms  as 
►  ,  '  _  z.  ,  ,  k(a  +  ib) 


£  +  in  =  (e  +  if)  + 


1  +  k(c  +  id) 


where 


f3  =  4  +  in 

p 


and  a,  b,  c,  d,  e,  and  f  are  all  real.  The  imaginary  part 
of  Equation  (8)  can  be  rearranged  to  the  equation  of  a 
circle,  being 


%  ^  b  \  1 2  .  T_  ic  a  \l2  _  a2  +  b2 

}  ~  le  2d  /  J  [n  "lf  '  2d  )  J  ”  4d 


It  is  important  to  note  here  that  the  radius  of  the  circular 
response,  as  well  as  the  complex  coordinates  of  the  center 
of  the  circle,  may  be  obtained  from  these  original  complex 


m 


receptance  matrix  coupling  terms  of  Equation  (10).  No 
iterations  in  k  are  necessary;  in  fact,  k  drops  out  of  the 
solution  when  expressed  in  this  form.  Furthermore,  it  is 
clearly  seen  that  these  circular  response  properties  are 
purely  a  function  of  the  imaginary  terms  (i.e.,  damping). 

Several  papers  have  been  published  (References  2,  3,  and 
1)  to  both  mathematically  and  experimentally  substantiate  this 
circular  response  region,  and  to  employ  these  circular  prop¬ 
erties  as  a  possible  optimi2ation  technique  for  reducing 
structural  vibration.  In  a  more  general  sense,  Balmford5 
describes  the  Vincent  circle  as  a  method  for  providing  a 
reduction  of  vibration  in  a  local  area  by  adjustment  of 
airframe  modes  such  that  modal  cancellation  will  take  place 
to  reduce  the  response  for  one  frequency  and  one  loading 
condition. 

STRAIN  ENERGY  METHOD 


Strain  energy  is  an  expression  for  the  potential  energy  of  a 
structural  element  and  is  most  commonly  expressed  in  matrix 
form  for  the  static  condition  as 

SE  =  1/2  (6TKe  6)  (11) 

For  the  dynamic  condition,  structural  elements  possessing  the 
highest  strain  energies  in  a  given  mode  of  vibration  have  been 
shown  to  be  the  best  candidates  for  modification  to  reduce 
overall  structural  dynamic  amplification. 6  Two  approaches  are 
generally  considered,  as  described  in  the  following  sections. 

Modal  Strain  Energy  Approach 

Tnis  approach  seeks  to  reduce  vibration  by  detuning  a  struc¬ 
ture  natural  frequency  to  a  better  position  with  respect  to 
the  forcing  frequency.  Modal  strain  energy  is  calculated 
for  each  structural  element  using  the  mode  shape  (eigenvector) 
for  the  natural  frequency  to  be  modified 

SE  =  1/2  (0TKe  0)  (12) 

Elements  with  the  highest  strain  energies  indicate  t:.e  optimal 
structural  elements  to  change  in  order  to  shift  the  natural 
frequency  a way  from  the  forcing  frequency,  thereby  reducing 
dynamic  amplification. 
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Forced  Response  Strain  Energy  Approach 

This  is  an  extension  of  the  modal  approach  that  takes  into 
account  the  response  participation  of  all  modes  for  a  parti¬ 
cular  load  application  and  excitation  frequency,  including 
the  effects  of  structural  damping  and  combined  forcing  func¬ 
tion  phase  relationships.  Therefore,  the  resultant  strain 
energy  of  each  structural  element  will  be  a  periodic  function 
in  time.  Sciatra6  has  done  considerable  work  in  this  field 
and  has  arrived  at  an  expression  for  the  maximum  damped  forced 
response  element  strain  energy  within  the  period 


4(SE)n,a*  =  6l  Ke6R  +  6>e6I 


+  [<6IKa6I  '  6RKe6R>2 


T  T  2  1/2 

+  <6RKe6I  +  6lW  3  7 


(13) 


Elements  with  the  highest  strain  energies  indicate  the  ele¬ 
ments  that  are  most  responsible  for  the  structural  dynamic 
amplification.  Because  the  forced  response  strain  energy  ap¬ 
proach  gives  direct  information  at  the  specific  excitation 
frequency  of  interest,  it  was  selected  for  comparison  to  the 
Vincent  circle  method. 


Figure  2  illustrates  an  application  of  the  strain  energy 
technique.  First,  the  modal  strain  energy  approach  was  used 
to  selectively  stiffen  the  structure  in  order  to  move  a 
fuselage  natural  frequency  to  a  higher  position  above  the 
excitation  frequency,  thereby  reducing  the  dynamic  amplifica¬ 
tion.  A  damped  forced  response  analysis  of  the  modified 
structure  was  then  compared  to  the  original  structure,  and 
it  was  confirmed  that  an  overall  vibration  reduction  had 
resulted  in  the  fuselage  area  of  interest. 
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Natural  Frequency 

9 

8.442 

8.494 

Changed  (Mode  No.  10, 

10 

12.03 

12. 74-. - 

4tn  Elastic  Mode)— 

11 

12.52 

12.79 

Mode  Amplification 

12 

14.48 

14.82 

Factor  Down  From 

13 

15.69 

16.02 

9.1  to  5.0 

14 

16.95 

17.22 

15 

18.94 

19.06 

A  damped  forced  response  analysis  comparison  of  the 
original  and  the  modified  fuselage  resulted  in  a 
vibration  reduction  for  99  out  of  123  degrees-of-f reedom. 


Figure  2.  Boeing-Vertol  CH-47A  vibration  reduction  by 
structural  modification  (Reference  6). 
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PROGRAM  OBJECTIVES 


The  purpose  of  this  contracted  study  was  to  evaluate  the 
practical  applicability  of  the  Vincent  circle  method,  as 
compared  to  the  forced  response  strain  energy  method,  for 
reducing  vibration  through  realistic  structural  stiffness 
optimization.  The  analytical  study  was  performed  on  the 
Model  AH-1G  helicopter,  and  the  effort  was  limited  to 
application  of  a  single  steady-state  main  rotor  2/rev 
(10.8  Hz)  vertical  excitation  force  with  a  goal  of  mini¬ 
mizing  the  pilot's  seat  vertical  response. 

The  Vincent  circle  method  was  further  investigated  for 
usefulness  as  applied  to  evaluating  mass  changes,  damping, 
and  ability  to  optimize  dynamic  absorber  locations. 
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ANALYTICAL  APPROACH 


METHOD  OF  ANALYSIS 


For  evaluating  the  Vincent  circle  method,  the  receptance 
matrix  technique  previously  described  offers  an  attractive 
and  direct  means  to  obtain  the  response  circle  radius  and 
center-of-location  properties  based  on  a  single  stiffness 
(or  mass)  parameter.  However,  for  a  feasibility  study  of 
this  nature,  several  important  disadvantages  also  exist  in 
that  the  receptance  matrix  technique  is  not  easily  adapted 
for  evaluation  of  the  following: 

The  sensitivity  of  the  circular  response  region  versus 
element  stiffness  (or  mass)  values. 

-  A  coupled  stiffness  matrix  between  two  points.  In 
dealing  with  real  structure,  rarely  are  we  concerned 
with  a  single  uncoupled  parameter  like  a  simple  spring, 
but  rather  structural  elements  behaving  more  like  beams 
that  have  coupled  stiffness  matrices.  For  example,  a 
single  variation  in  the  area  moment  of  inertia  (I)  of  a 
beam  affects  both  the  shear  translational  and  rotational 
degrees-of-freedom  between  the  two  ends  of  the  beam. 

Individual  element  damping  variations  while  the  damping 
for  the  remainder  of  the  structure  remains  uniform. 

Tuning  sensitivity  of  a  dynamic  absorber  at  different 
locations  on  the  structure  (i.e.,  effect  of  backup 
structure  at  absorber  location). 

Since  the  purpose  of  this  contracted  effort  was  not  to  evalu¬ 
ate  the  merits  of  the  receptance  matrix  technique,  and  in 
order  to  use  the  same  analytical  tool  for  evaluating  both  the 
Vincent  circle  method  and  the  forced  response  strain  energy 
method,  a  straightforward  linear  NAS TRAN'  analysis  was  se¬ 
lected  for  performing  the  study.  Also,  due  to  a  unique 
situation  at  Bell  Helicopter  Textron  wherein  many  NASTRAN 
computer  runs  can  be  accomplished  quickly  and  at  a  low  cost, 
the  additional  time  and  cost  that  would  have  been  required 
to  develop  a  special  purpose  receptance  matrix  manipulation 
computer  program  was  avoided. 


7THE  NASTRAN  USER'S  MANUAL,  NASA  SP-222(03)  National  Aero¬ 
nautics  and  Space  Administration,  Washington,  D.  C.,  July 
1976. 
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As  applied  to  the  Vincent  circle  method,  this  straightforward 
NASTRAN  approach  required  the  evaluation  of  several  iterations 
in  the  particular  parameter  under  investigation  to  develop 
the  circular  response  region,  from  which  the  response  circle 
radius  and  center-of-location  properties  were  subsequently 
determined.  However,  no  restrictions  were  necessary  on  what 
parameters  could  be  selected,  and  more  in-depth  information 
was  provided  as  to  the  characteristics  of  the  circular  re¬ 
sponse  region  than  if  the  receptance  matrix  technique  had  been 
employed. 

The  version  of  NASTRAN  used  at  the  time  of  this  study  was 
Level  16.0.5.  This  particular  level  of  NASTRAN  does  support 
strain  energy  calculations,  but  only  for  static  analysis. 
Although  desirable,  it  was  found  that  the  NASTRAN  internal 
modifications  necessary  to  incorporate  the  complex  damped 
forced  response  strain  energy  calculations  (Equation  13)  were 
beyond  the  scope  of  this  study.  However,  by  using  the  Direct 
Matrix  Abstraction  Program  (DMAP)  capability  in  NASTRAN,  the 
static  strain  energy  module  was  incorporated  into  the  dynamics 
frequency  response  analysis  with  the  following  limitations: 

All  applied  oscillatory  loads  must  be  input  at  zero 
degrees  phase. 

-  No  damping  can  be  included,  thus  the  applied  load  and 
response  vector  have  only  real  components. 

These  limitations  were  required  because  the  statics  format 
of  the  NASTRAN  strain  energy  module  is  not  compatible  for 
processing  complex  numbers  (i.e.,  only  the  real  part  of  a 
complex  number  is  considered).  Therefore,  the  forced  re¬ 
sponse  strain  energy  expression  for  a  structure  element  re¬ 
sponding  to  a  single  resultant  applied  load  with  zero  damp¬ 
ing  has  now  been  reduced  to 

3E  =  1/2  («Jk66r)  (14) 


This  results  in  a  strain  energy  distribution  that  is  a  close 
approximation  to  the  strain  energy  distribution  obtained  when 
a  small  amount  of  uniform  structural  damping  is  included.  For 
the  purposes  of  this  study,  this  undamped  strain  energy  dis¬ 
tribution  was  acceptable  for  comparison  to  the  Vincent  circle 
results,  wherein  a  small  amount  of  damping  must  be  included  as 
previously  explained. 
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MATHEMATICAL  MODELS 


Initial  stiffness  parameter  investigations  were  conducted 
using  a  simple  elastic-line  NASTRAN  model  of  the  AH-1G  heli¬ 
copter  (Figures  3  and  4).  This  elastic-line  model  was  originally 
developed  as  an  AH-1J  (Reference  8)  and  was  later  modified 
into  the  AH-1G  configuration  (Reference  9).  The  structural 
optimization  method  showing  the  most  potential  for  reducing 
vibration  through  realistic  structural  stiffness  changes  was 
then  applied  to  minimize  pilot's  seat  vertical  response  using 
a  large  complex  built-up  NASTRAN  model  of  the  AH-1G  (Figure 
5).8 *  10  The  elastic-line  model  was  also  used  to  further  evalu¬ 
ate  the  Vincent  circle  method  as  applied  to  mass,  damping,  and 
dynamic  absorber  parameters . 

For  use  in  this  study,  both  the  elastic-line  and  the  built- 
up  models  of  References  9  and  10,  respectively,  were  updated 
to  a  basic  mission-clean  wing  configuration.  To  show  the 
similarity  of  the  two  models,  a  comparison  of  weight  data  and 
natural  frequency  placements  is  shown  in  Table  1.  From  the 
natural  frequency  comparisons,  in  addition  to  the  obvious 
absence  of  landing  gear  skid  modes,  it  is  seen  that  the 
elastic-line  model  has  no  fuselage  torsional  modes  in  the 
zero  to  thirty  hertz  frequency  range  of  interest.  This  is 
due  to  inadequate  fuselage  torsional  mass  inertia  representa¬ 
tion  in  the  elastic-line  model  which,  due  to  the  torsional 
coupling,  is  also  partially  responsible  for  the  differences 
in  pylon  roll  and  fuselage  lateral  bending  modes.  Note, 
however,  the  close  agreement  between  the  two  models  for 
the  pylon  pitch  and  the  first  and  second  fuselage  vertical 
bending  modes  that  are  of  primary  importance  in  this  study. 

Both  the  elastic-line  and  built-up  models  have  the  same 
elastic-line  representations  for  the  helicopter  tailboom 
structure . 


8Cronkhite,  J.  D. ,  and  Wilson,  W.  F.,  DYNAMIC  ANALYSIS  OF 
TWO-PER-REV  VIBRATIONS  IN  THE  MODEL  AH-1J  HELICOPTER  -  PIP 
Task  No.  AH-8-123,  Bell  Helicopter  Textron  Report  299-100- 
021,  Fort  Worth,  Texas,  4  February  1972. 

aCronkhite ,  J.  D. ,  XM-97  (20MM)  WEAPON  ON  THE  AH-1G  -  PRELI¬ 
MINARY  DYNAMIC  ANALYSIS,  Bell  Helicopter  Textron  Inter¬ 
office  Memo  81 : JDC : mb -054 ,  Fort  Worth,  Texas,  29  May  1973. 

10Cronkhite,  J.  D.,  Berry,  V.  L.,  and  Brunken,  J.  E.,  A 
NASTRAN  VIBRATION  MODEL  OF  THE  AH-lG  HELICOPTER  AIRFRAME, 
U.S.  Army  Armament  Command  Report  No.  R-TR-64-45,  Re¬ 
search  Directorate,  Gen.  Thomas  J.  Rodman  Laboratory, 

Rock  Island  Arsenal,  Rock  Island,  Illinois,  June  1974. 
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Figure  3.  AH-1G  elastic-line  math  model. 


Figure  4. 


Elastic-line  model  main  rotor  pylon 
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Figure  5.  AH- 1 G  built-up  math  model. 


TABLE  1.  AH- 1 G  MODEL  COMPARISON 


Basic  Mission-clean  Wing 

3-D 

Built-Up 

Model 

Elastic 

Line 

Model 

Weight  Data 

Gross  weight  (lb) 

8394  . 

8394. 

{ Sta  (in.) 

192.9 

193.1 

Center-of-gravity  <  BL  (in.) 

0.0 

0.1 

(wL  (in.) 

70.6 

68.7 

/ Roll  (lb-in . 2 ) 

1.139  x  107 

1.139  x  107 

Inertias  < Pitch  (lb-in.2) 

6.193  x  107 

6.200  x  107 

(Yaw  (lb-in.2) 

5.196  x  107 

5.215  x  107 

Natural  Frequency  Data 

Mode  Description 

(Hz) 

(Hz) 

M/R  pylon  pitch 

3.03 

3.02 

M/R  pylon  roll 

3.90 

4.24 

1st  Fuselage  lateral  bending 

7.14 

6.80 

1st  Fuselage  vertical  bending 

7.94 

7.93 

1st  Skid  mode 

14.57 

- 

1st  Fuselage  torsion 

15.66 

- 

2nd  Fuselage  lateral  bending 

17.49 

16.70 

2nd  Fuselage  vertical  bending 

17.49 

17.86 

2nd  Skid  mode 

18.76 

- 

3rd  Skid  mode 

19.84 

— 

2nd  Fuselage  torsion 

21.49 

- 

4th  Skid  mode 

23.43 

- 

M/R  mast  lateral  bending 

25.28 

24. 79 

M/R  mast  F/A  bending 

24 .97 

25.80 

5th  Skid  mode 

25.75 

- 

3rd  Fuselage  vertical  bending 

26.96 

29.47 

6th  skid  mode 

29.04 

— 

The  AH-lG  elastic-line  NASTRAN  model  (including  the  DMAP  ALTER 
procedure  developed  for  obtaining  undamped  forced  response 
element  strain  energy  output)  is  included  as  Appendix  A. 

PRACTICAL  CRITERION  FOR  STRUCTURAL  STIFFNESS  CHANGES 


A  practical  criterion  for  stiffness  change  limitations  was 
developed  by  considering  the  classical  single  degree-of- 
freedom  spring/mass/damper  dynamic  system  that  has  a  re¬ 
sonance  dynamic  amplification  factor  of  25  for  2  percent 
damping.  This  amplification  factor  can  be  reduced  to  1.1 
by  either  reducing  the  stiffness  by  a  factor  of  0.5,  or 
by  increasing  the  stiffness  by  a  factor  of  10,  assuming  the 
mass  and  damping  parameters  to  remain  constant.  This  0.5 
to  10  stiffness  factor  range  seems  to  be  within  the  limits 
of  practical  considerations  even  though,  realistically 
speaking,  structure  designed  for  strength  would  not  have 
its  stiffness  reduced  by  one-half  or  increased  ten-fold 
without  undergoing  a  considerable  redesign  effort  (i.e., 
geometric  shape  factor,  type  materials,  added  weight 
penalties,  etc. ) 


Based  on  this  practical  stiffness  change  criterion,  the 
forced  response  strain  energy  and  Vincent  circle  results 
were  compared  to  determine  which  was  the  most  promising 
method  for  vibration  reduction  through  realistic  structural 
stiffness  optimization. 
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STIFFNESS  PARAMETER  INVESTIGATIONS 


ELASTIC-LINE  MODEL  ANALYSIS 


The  forced  response  strain  energy,  method  was  used  to  deter¬ 
mine  which  structural  members  had  the  highest  forced  re¬ 
sponse  strain  energies.  Using  stiffness  parameter  varia¬ 
tions,  these  members  were  then  analyzed  using  the  Vincent 
circle  method. 

Forced  Response  Strain  Energy  Results 

Figure  6  is  a  NASTRAN-generated  plot  showing  a  side  view  of 
the  elastic-line  model  forced  response  vertical  deformation 
mode  shape.  The  resultant  forced  response  strain  energy 
distribution  is  shown  in  Figure  7.  The  element  identifica¬ 
tion  numbers  shown  are  related  to  the  two-digit  GRID  num¬ 
bers  at  each  end  of  the  element,  as  identified  in  Figures  3 
and  4.  The  strain  energy  shown  for  each  element  is  the 
total  due  to  the  resultant  six  degrees-of-freedom  deflec¬ 
tions  at  each  end  of  the  element  and  the  element  stiffness 
matrix  as  calculated  by  Equation  (14).  Since  the  elastic¬ 
line  model  is  essentially  symmetric  and  the  applied  load 
at  the  hub  is  in  the  vertical  direction,  only  deflections 
in  the  vertical  plane  are  contributing  significantly  to  the 
strain  energy  calculation.  Note  how  the  strain  energy  dis¬ 
tribution  compares  to  what  might  be  expected  from  examina¬ 
tion  of  the  elastic  deformation  mode  shape  of  Figure  6. 

The  elastic-line  model  NASTRAN  forced  response  strain  energy 
output  data  are  presented  in  Appendix  B. 

Vincent  Circle  Results 


As  previously  shown,  this  circular  response  property  is  based 
on  variations  involving  single  parameters  only.  However, 
in  this  study  the  Vincent  circle  method  was  evaluated  with 
respect  to  realistic  element  property  changes  such  as  AE 
axial  stiffness  and  El  bending  stiffness  parameters,  wherever 
applicable,  so  that  the  influence  of  all  stiffness  coupling 
terms  would  be  included. 

In  all  cases,  a  single  1000-pound  2/rev  main  rotor  vertical 
excitation  force  was  applied  and  the  vertical  response  of 
the  pilot's  seat  was  calculated.  Unless  otherwise  stated, 
all  Vincent  circle  calculations  were  performed  assuming  2- 
percent  uniform  structural  damping. 
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model  2/rev  forced  response  mode  shape. 
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Figure  7.  Elastic-line  model  forced  response  strain  energy- 
distribution. 
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The  first  application  of  this  method  was  to  verify  the  cir¬ 
cular  response  property  using  the  elastic-line  model.  This 
was  accomplished  by  evaluation  of  the  AE  axial  stiffness 
parameter  for  the  upper  main  rotor  mast  segment,  as  shown  in 
Figure  8.  Of  particular  significance  are  the  zero  response 
at  zero  stiffness  and  the  sensitivity  of  response  to  stiffness 
variations.  The  zero  response  at  zero  stiffness  occurs  since 
the  main  rotor  mast  provides  the  only  load  path.  The  baseline 
stiffness  value  of  the  element  indicates  a  low  initial  amplifi¬ 
cation  factor,  and  the  response  is  seen  to  be  relatively 
insensitive  to  stiffness  variations  except  for  drastic  changes 
over  a  very  small  range  of  values  where,  in  reality,  the  element 
stiffness  is  the  single  parameter  responsible  for  the  resonance 
condition.  Since  the  applied  load  was  directed  straight  down 
the  mast  axis,  only  the  mast  AE  stiffness  parameter  was 
evaluated. 

With  the  circular  response  property  thus  verified,  the  re¬ 
mainder  of  the  thirteen  elements  showing  significant  forced 
response  strain  energies  were  evaluated  for  Vincent  circle 
properties  of  circle  diameter  and  center  location.  In  each 
case,  all  parameters  were  maintained  at  their  original  base¬ 
line  values  except  for  the  particular  stiffness  parameter 
under  investigation.  Aft  fuselage  and  tailboom  elements 
were  evaluated  separately  for  both  AE  and  vertical  El  stiff¬ 
ness  parameter  variations.  It  was  found  that  element  AE 
parameter  changes  produced  a  response  circle  which  included 
the  baseline  AE  stiffness  response  point,  as  previously  shown 
for  element  5455  in  Figure  8.  Changes  in  the  element  El 
parameter  were  found  to  produce  another  response  circle,  but 
one  that  was  offset  from  the  baseline  El  stiffness  response 
point.  Further  investigations  showed  that  this  offset  of 
the  El  response  circle  was  due  to  the  I  associated  stiffness 
coupling  between  beam  transverse  deflection  and  beam  bending 
(slope  change)  degrees-of-freedom.  This  stiffness  coupling 
offset  effect  was  found  to  be  most  pronounced  for  element 
1415  (Figure  9).  The  pylon  lift  link,  element  5112,  and  the 
pylon  5th  mount,  element  4550,  are  analogous  to  axially  loaded 
pin-ended  rods  so  that  only  AE  stiffness  parameters  were 
evaluated. 

The  response  circle  diameter  indicates  the  maximum  response 
change  possible  due  to  a  particular  element  stiffness  param¬ 
eter.  A  normalized  circle  diameter  distribution  is  shown 
in  Figure  10.  Figure  11  shows  the  maximum  reduction  in  re¬ 
sponse  due  to  each  element  stiffness  parameter  as  determined 
from  its  circle  diameter  and  center  location.  In  the  case 
of  elements  evaluated  for  both  AE  and  El  stiffness  parameters, 
the  parameter  providing  the  maximum  reduction  in  response  was 
selected.  Note  that  the  trend  of  element  effectiveness  estab¬ 
lished  in  Figure  10  is  similar  to  that  depicted  in  Figure  11. 
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Simple  elastic  line  model  with  2%  uniform  structural  dampin 


Figure  9.  Vincent  circle  plots  for  aft  fuselage  segment. 


Normalized  to  largest  circle  diameter 
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In  other  words,  in  this  application  the  elements  with  the 
largest  response  circle  diameters  are  also  capable  of  achiev¬ 
ing  the  minimum  response;  however,  this  is  not  to  be  con¬ 
strued  as  true  in  all  cases. 

The  most  important  result  to  be  pointed  out  is  that  the 
Vincent  circle  properties  of  circle  diameter  and  center-of- 
location  in  the  complex  plane  clearly  can  be  used  to  deter¬ 
mine  the  minimum  response  attainable  for  a  particular  stiff¬ 
ness  parameter;  however,  no  information  is  provided  as  to 
what  stiffness  value  is  required  to  achieve  this  minimum 
response.  Table  2  shows  the  drastic  and  unrealistic  stiff¬ 
ness  values  that  were  required  to  produce  both  the  minimum 
and  the  maximum  pilot's  seat  vertical  response  points  on 
the  response  circle  defined  by  each  element  stiffness  param¬ 
eter.  The  two  cases  where  negative  stiffnesses  are  indi¬ 
cated  occurred  because  several  elements  provided  parallel 
load  paths  so  that  a  negative  stiffness  in  one  element  was 
required  to  sufficiently  reduce  the  combined  elements  posi¬ 
tive  stiffness  resultant.  In  most  cases,  there  is  very 
little  difference  between  the  element  stiffness  values  at 
minimum  response  and  at  maximum  response  (with  exceptions 
for  elements  1415,  1516,  and  1617  due  to  their  El  response 
circle  offsets  (see  Figure  9)).  Figure  12  further  illus¬ 
trates  the  stiffness  nonlinearity  around  the  response  circle 
by  depicting  the  typical  stiffnesses  generally  found  to  des¬ 
cribe  the  majority  of  the  circular  response  region.  This 
inherent  property  can  be  somewhat  rationalized  by  considering 
a  single  degree-of- freedom  model  where  the  dynamic  amplifica¬ 
tion  is  greatest  for  the  resonance  condition  and  decreases 
rapidly  as  the  natural  frequency  is  shifted  either  higher  or 
lower  than  the  excitation  frequency. 

Strain  Energy/Vincent  Circle  Comparison 

Figure  13  compares  the  undamped  forced  response  strain 
energy  distribution  to  the  2-percent  damped  Vincent  circle 
diameter  distribution  for  the  same  elements  of  the  elastic¬ 
line  model.  Here  again,  where  multiple  stiffness  parameters 
were  evaluated  for  the  same  element,  the  stiffness  parameter 
producing  the  largest  circle  diameter  was  chosen.  It  is 
seen  that  the  Vincent  circle  method  does  not  give  the  same 
distribution  picture  as  the  forced  response  strain  energy 
method.  The  forced  response  strain  energy  method  points 
to  the  tailboom  as  the  area  most  responsible  for  dynamic 
amplification  due  to  elastic  deformation,  while  the  Vincent 
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TABLE  2.  ELASTIC-LINE  MODEL  ELEMENT  STIFFNESS  VALUES 

FOR  MINIMUM  AND  MAXIMUM  PILOT'S  SEAT  RESPONSE 


Stiffness 

%  of  Baseline  Stiffness  | 

Element 

Parameter 

Minimum 

Maximum 

No . 

Type 

Response 

Response 

1415 

AE 

.  021% 

.  024% 

El 

100% 

.  071% 

1516 

AE 

.  011% 

.  013% 

El 

100% 

.  027% 

1617 

AE 

.  190% 

.  213% 

El 

100% 

.  015% 

1718 

AE 

.  181% 

.  206% 

El 

.  020% 

.  02 1% 

1819 

AE 

.  190% 

.  210% 

El 

.  024% 

.  028% 

1920 

AE 

.  168% 

.  188% 

El 

.  027% 

.  031% 

2021 

AE 

.152% 

.  168% 

El 

.  035% 

.  041% 

2122 

AE 

.123% 

.139% 

El 

.  036% 

.  043% 

2223 

AE 

.115% 

.  132% 

El 

.  056% 

.  071% 

2324 

AE 

.104% 

.  118% 

El 

.  109% 

.190% 

51]  2 

AE 

-.775% 

-.570% 

5455 

AE 

0% 

.  387% 

4550 

AE 

-167% 

-147% 

K\X\\  ,'l  Strain  energy  distribution  normalized  to  100%  total 

!'; . I  Largest  Vincent  circle  diameter  per  element  normalized 

to  100%  total 
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circle  method  points  to  the  pylon  as  the  area  having  the  most 
potential  for  reducing  vibration  at  the  pilot's  seat. 

The  pylon  elements  indicated  by  the  Vincent  circle  method  to 
have  the  greatest  potential  for  reducing  vibration  were  further 
evaluated.  It  was  found,  as  might  be  strongly  suspected  from 
the  previous  discussion  of  results,  that  stiffness  variations 
for  these  pylon  elements  within  the  established  confines  of  the 
0.5  to  10  stiffness  factor  range,  even  if  done  collectively, 
yielded  virtually  no  substantial  change  in  response  at  the 
pilot's  seat.  However,  stiffness  changes  of  this  same  magni¬ 
tude  applied  to  those  tailboom  elements  initially  undergoing 
substantial  elastic  deformation,  as  indicated  by  the  forced 
response  strain  energy  method,  were  found  to  produce  signi¬ 
ficant  changes  in  response  at  the  pilot's  seat.  Therefore, 
the  forced  response  strain  energy  method  was  determined  to  be 
better  adapted  for  vibration  reduction  through  structural 
stiffness  optimization  by  indicating  which  structural  elements 
are  responsible  for  the  initial  dynamic  amplification;  it  was 
also  determined  that  realistic  structural  stiffness  changes  in 
these  elements  can  efficiently  alter  the  vibration  character¬ 
istics  of  the  structure. 

BUILT-UP  MODEL  VIBRATION  OPTIMIZATION 

At  this  point  in  the  investigations,  the  selected  forced 
response  strain  energy  method  was  applied  to  a  large  complex 
built-up  NASTRAN  model  (Figure  5). 

A  NASTRAN-generated  plot  of  a  side  view  of  the  built-up  model 
forced  response  vertical  deformation  mode  shape  is  shown  in 
Figure  14.  The  resultant  undamped  forced  response  strain 
energy  distribution  as  compared  to  the  distribution  for  the 
elastic-line  model  is  shown  in  Figure  15.  The  elastic-line 
model  element  identification  numbers  were  retained  on  the 
figure  for  reference.  Note  the  similarity  of  the  two  distri¬ 
butions.  In  the  case  of  the  aft  fuselage  and  main  rotor  mast 
elements,  the  single  element  representation  of  the  elastic¬ 
line  model  has  been  replaced  in  the  built-up  model  by  more 
detailed  modeling  involving  many  elements  more  uniformly  shar¬ 
ing  the  load,  and  each  accounts  for  a  smaller  percentage  of 
the  total  strain  energy. 

Previous  analyses6  have  verified  that  the  most  dynami  .ally 
efficient  structure  for  a  given  mode  of  vibration  is  one  with 
a  uniform  energy  distribution.  Therefore,  the  most  efficient 
structural  stiffness  modifications  incorporated  for  pilot's 
seat  vibration  reduction  should  also  result  in  a  more  uniform 
strain  energy  distribution. 
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Figure  15.  Comparison  of  forced  response  strain  energy 
distributions . 
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From  the  baseline  forced  response  strain  energy  distribution, 
the  tailboom  elements  are  seen  to  possess  the  most  strain 
energy  and  were  selected  as  the  best  candidates  for  stiffness 
modifications.  The  first  step  was  to  evaluate  both  extremes 
of  the  selected  practical  stiffness  criterion  range.  For  the 
stiffer  case,  the  tailboom  element  with  the  most  strain  energy 
was  made  ten  times  stiffer,  and  the  other  tailboom  elements 
were  made  stiffer  according  to  their  baseline  strain  energy 
ratio.  For  the  softer  case,  the  tailboom  element  with  the 
least  strain  energy  was  made  one-half  as  stiff,  and  the  other 
tailboom  elements  were  reduced  in  stiffness  according  to  their 
baseline  strain  energy  ratio  (an  additional  criterion  had  to 
be  imposed  here  -  the  resultant  element  stiffness  could  not 
be  increased  because  of  this  strain  energy  ratio  scheme  used). 
The  strain  energy  ratio  scheme  is  shown  in  further  detail  in 
Table  3,  along  with  the  accompanying  results.  In  both  of 
these  cases  the  resultant  strain  energies  are  more  uniformly 
distributed.  In  both  cases,  the  pilot's  seat  response 
increased  in  magnitude,  but  for  the  stiffer  case  the  response 
changed  phase.  This  result  indicates  that  a  proper  stiffness 
increase,  being  within  the  selected  practical  stiffness 
criterion  range,  would  result  in  an  absolute  zero  response  at 
the  pilot's  seat.  Two  additional  stiffness  iterations  were 
required  to  achieve  the  desired  results.  In  each  case  the 
selected  stiffness  ratio  was  applied  to  tailboom  element 
number  4014220  (largest  baseline  strain  energy)  and  the  other 
tailboom  elements  were  made  stiffer  according  to  their  base¬ 
line  strain  energy  ratio.  As  shown  in  Table  3,  it  was  found 
that  a  maximum  stiffness  increase  of  375  percent  resulted  in 
near  zero  response  at  the  pilot's  seat  and  reduced  the  total 
strain  energy  in  the  tailboom  by  46  percent.  This  also 
resulted  in  a  slight  increase  in  the  strain  energies  of  the 
two  tailboom  elements  nearest  the  fuselage.  Additional  selec¬ 
tive  stiffening  of  the  tailboom  could,  no  doubt,  produce  the 
same  zero  response  results  with  a  more  uniform  strain  energy 
distribution;  however,  as  it  was  not  the  purpose  of  this  study 
to  develop  a  general  optimization  procedure,  further  investi¬ 
gations  in  this  area  were  not  pursued. 

A  better  interpretation  of  these  results  may  be  gained  from 
examination  of  the  simplified  forced  response  mode  shapes 
shown  in  Figure  16.  Here  it  can  easily  be  seen  how  the 
deformed  mode  shape  was  altered  through  the  stiffness  modi¬ 
fications  to  place  the  pilot's  seat  at  a  node  point. 

Although  in  the  foregoing  analysis  it  was  possible  to  achieve 
zero  dynamic  response  at  a  particular  point  under  a  unique  set 
of  circumstances,  in  no  way  should  these  results  be  inter¬ 
preted  as  implying  that  zero  response  is  always  attainable 


lqure  J6.  Simplified  built-up  model  forced  response  mode 
shapes  showing  vibration  optimization  process. 


through  any  type  of  structural  modification.  As  previously 
explained  in  the  introduction,  many  factors  need  to  be  con¬ 
sidered.  Also,  for  the  purposes  of  this  study,  no  change  in 
mass  was  considered  for  any  of  the  stiffness  changes  evaluated. 
One  of  the  next  most  logical  steps  necessary  for  any  realistic 
vibration  optimization  program  would  be  to  evaluate  the  com¬ 
bined  effects  of  both  stiffness  and  the  associated  weight 
change.  The  end  product  is,  of  course,  a  minimum  vibration 
environment  obtained  through  optimization  of  a  structure 
within  a  given  set  of  practical  criteria. 


45 


FURTHER  EVALUATION  OF  THE  VINCENT  CIRCLE 


The  Vincent  circle  phenomenon  was  further  evaluated  for  mass 
tuning,  damping,  and  dynamic  absorber  parameters  using  the 
elastic-line  model  and  a  simple  add-on  single  degree-of- 
freedom  spring-mass-damper  system. 

MASS  PARAMETER  RESULTS 

Element  mass  parameter  changes  were  found  to  describe  the  same 
circle  that  stiffness  changes  produce,  only  the  increasing 
mass  traced  out  the  circle  in  the  opposite  direction.  Data 
points  on  the  example  shown  in  Figure  17  have  been  somewhat 
exaggerated  for  clarity. 

DAMPING  PARAMETER  RESULTS 

Damping  is  usually  not  considered  as  a  practical  method  of 
effectively  altering  airframe  structural  vibration.  As 
shown  in  Figure  18,  both  the  maximum  amplification  and  the 
potential  of  the  element  for  vibration  control  are  greatly 
reduced  by  small  increases  in  element  damping. 

DYNAMIC  ABSORBER  PARAMETER  RESULTS 

The  adopted  procedure  for  investigating  dynamic  absorber 
parameters  was  to  evaluate  the  effectiveness  of  one  absorber 
for  attenuating  2/rev  vertical  vibration  at  the  pilot's 
seat  when  placed  at  different  locations  on  the  structure; 
and  then,  to  compare  these  results  to  the  maximum  reductions 
in  pilot's  seat  response  obtained  from  element  stiffness 
Vincent  circle  properties  at  these  same  locations  (as  des¬ 
cribed  in  Figure  11 ) .  The  adopted  approach  does  not  imply 
that  any  such  relation  between  Vincent  circle  predictions  and 
optimum  dynamic  absorber  locations  exists,  or  has  been 
reported  to  exist. 

The  dynamic  absorber  energy  absorption  potential  was  kept 
constant  by  maintaining  a  25-pound  absorber  mass  and  a  0.2% 
absorber  damping  coefficient.  This  required  several  absorber 
spring  stiffness  iterations  to  achieve  the  proper  tuning  of 
the  dynamic  absorber  at  each  location  due  to  the  different 
boundary  conditions  (i.e.,  backup  structure  stiffness)  which 
exists  throughout  the  structure.  In  each  case  the  dynamic 
absorber  was  positioned  in  the  vertical  direction  for  re¬ 
ducing  pilot's  seat  vertical  vibration.  For  this  study 
seven  dynamic  absorber  locations  were  evaluated  and  the 
results  are  shown  in  Figure  19  as  compared  to  the  maximum 
response  reduction  distribution  obtained  from  the  element 
stiffness  Vincent  circle  properties. 
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Maximum  reduction  in  response  based  on 
element  stiffness  Vincent  circle  properties 


12/(T)  8  9/0  1516/0  1920/O  2324/(0)  5112/0  5455/© 

Element/(^rid)  identification  numbers 


Nose  Pilot's  Fwd  Mid  Aft  Pylon  Upper 

seat  T/B  T/B  T/B  lift  M/R 

link  mast 


Figure  19.  Comparison  of  Vincent  circle  properties  to 
dynamic  absorber  location  effectiveness  for 
reducing  pilot's  seat  vibration. 
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Examination  of  the  dynamic  absorber  location  effectiveness 
shows  that  the  most  efficient  location  for  reducing  pilot's 
seat  vibration  is  at  the  pilot's  seat;  but,  it  also  shows  that 
a  remote  absorber  located  either  at  the  top  of  the  main  rotor 
mast  (grid  55),  at  the  pylon  lift  link  (grid  12),  or  at  the 
middle  of  the  tailboom  (grid  20),  would  also  be  significantly 
effective  for  reducing  vibration  at  the  pilot's  seat. 

Use  of  the  Vincent  circle  properties  shows  that  both  elements  89 
and  5455  are  capable  of  reducing  the  pilot's  seat  vibration  by 
100%.  However,  these  results  are  inconclusive  because  any 
element  providing  the  only  load  path  is  ultimately  capable  of 
reducing  the  pilot's  seat  response  to  zero  if  its  stiffness  is 
reduced  to  zero.  Further  examination  of  the  comparisons  shown 
in  Figure  19  shows  no  definite  trends  to  indicate  that  the 
Vincent  circle  method  can  be  applied  to  predict  optimum  dynamic 
absorber  locations . 
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CONCLUSIONS 


The  Vincent  circle  property  of  a  circular  response  region  has 
been  examined  and  verified  by  .straightforward  linear  analysis. 
This  property  would  certainly  be  of  value  to  the  dynamicist 
for  evaluation  of  a  localized  portion  of  structure  known  to  be 
in,  or  very  near,  resonance  and  which  could  be  controlled  by 
changes  involving  only  a  few  parameters . 

However,  for  vibration  reduction  through  optimization  of  an 
airframe  structure  where  many  elements  are  involved,  the  forced 
response  strain  energy  method  was  found  to  be  more  suitable  than 
the  Vincent  circle  method  by  indicating  which  structural  elements 
are  most  responsible  for  the  elastic  dynamic  amplification. 
Practical  parameter  modifications  of  these  elements  will  most 
significantly  alter  the  airframe  structural  dynamic  charac¬ 
teristics  to  reduce  vibration.  For  the  more  detailed  analysis 
involving  complex  applied  loads,  phasing,  and  damping,  the 
damped  forced  response  strain  energy  method  has  the  potential 
benefit  of  rapid  evaluation  and  optimization  of  structural 
vibratory  response,  which  would  be  particularly  beneficial  in 
the  predesign  stages  of  an  aircraft. 
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CP  -  MUST  HE  A  I  — C  OL  UM  N 
TTfc*  FIRST  RUW=()..  A 
HP  -  MUST  BE  A  I  -COLUMN 
ALL  HU«fS=l.  ,  FXCFPT 


UENClF.S  COMPUTED 
S  MUST  OF  INPUT 

MY  NTR^Q-ROW  MATRIX. 

LL  UTHE^  IMTS-U 
MY  NIRc.vl-R»*A  matrix, 
fur  L  As  f  Kl  I  W  -  0. 

JPPUUm  F*  M  1*  !•***!* 
i*i*t*i*  inipunn  i  %  *  * 


Mi  rnm* 

mu  a  sir. 


0  00  I 
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GW  ID 

1 

39  .5  0 

0.  0 

50.0 

GR  ID 

2 

46  .0 

0.  0 

50 .0 

GW  ID 

3 

61  .2  5 

0.  0 

50  .0 

GW  ID 

4 

75.0 

0.  0 

50.0 

GR  ID 

5 

75.5 

0.  0 

29.0 

20  MM  GUN 

GW  ID 

6 

93  .0 

0.  0 

50 .0 

GW  ID 

7 

11  3.  5 

0.0 

5<J  .0 

GW  ID 

8 

13  1.5 

0.  0 

50 .0 

PILOT 

(Ji  ID 

y 

140.  0 

0.  0 

50  .0 

GW  ID 

10 

148.  5 

0.  0 

50.0 

GW  ID 

ti 

190.  0 

0.  0 

50  .0 

GW  ID 

12 

197.  0 

0.  0 

50.0 

GR  ID 

13 

21  3.  0 

0.  0 

50  .0 

GW  ID 

14 

25  0.  0 

O.  0 

50 .0 

GW  ID 

15 

27  7.  0 

0.0 

50.0 

GW  ID 

16 

299.  57 

0.  0 

49.59 

GW  ID 

17 

31  7.  72 

0.  0 

50 .7  2 

GW  ID 

1ft 

338.  6 1 

0.  0 

52.0  1 

GW  ID 

19 

35  9.  51 

0.  0 

53.3  I 

GW  ID 

20 

38  0.  42 

0.  0 

54  .6  1 

GW  ID 

21 

40  1.  JJ 

0*  0 

55  .9  0 
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22 

422.  24 

0.  0 

57.20 
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23 

44  J.  15 

0.  0 

58.5  0 
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24 

46  4.  10 

0.  0 

59.8  0 
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25 

48  0.  23 

0.  0 

75.4  0 
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2b 

48  8.  93 

0.  0 
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49  7.  77 
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0.  0 

109.  46 

GW  ID 

30 

52  0.  67 

0.  0 

11  ft.  27 

GW  ID 

31 

520.  67 

7.  20 

1 1  ft. 27 

GW  ID 

32 

52  0.  67 
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1000. 
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200. 
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1000. 
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o  r 

1 

1  000  . 
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2750. 

850. 
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1000. 

7300. 
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2000. 
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10  11 

1 
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12750. 
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1 
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14  71 0. 
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12  13 

1 

1  000  . 
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1  3  14 

1 

1  000  . 

1 3900. 

11670. 

2900. 

PH  AM 

14  16 

1 

1  ooo. 

I 21 00. 

I 0400. 

2650. 

PHAM 

1  6  IS 

1 

1  ooo  • 

0850  . 

9300  • 

2100. 

PHAM 

10  17 

1 

42. 7403 

5130. 3  9 

6244  .61 

3326.62 

PH  AM 

17  18 

1 

4  7.4632 

4606.89 

6300.05 

2307.82 

PRAM 

18  19 

1 

4 1 .4b92 

33  76  .8  0 

4650 .61 

1063.98 

PHAM 

l  920 

1 

38. 7424 

2797.46 

3323.30 

1422.92 

PH  AM 

2  021 

1 

3  7. 5l>06 

24 19.04 

230 1 .92 

1113.34 

PHAM 

2  1  22 

1 

36. 7154 

1919.77 

1782 . 38 

801.70 

PHAM 

22  23 

1 

34.7744 
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1 
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APPENDIX  B 

NASTRAN  FORCED  RESPONSE  STRAIN  ENERGY  OUTPUT  DATA 
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SYMBOLS 


A  Element  cross-sectional  area 

C  Viscous  damping  matrix 

E  Modulus  of  elasticity 

F  Applied  force 

G  Complex  receptance  matrix 

I  Element  area  moment  of  inertia 

K  Stiffness  matrix 

k  Scalar  spring 

M  Mass  matrix 

m  Scalar  mass 

S  Stiffness  ratio 

SE  Strain  energy 

t  Time 

x  Structural  displacement  response  vector 

Z  Pilot's  seat  vertical  displacement  response 

<5  Element  displacement  response  vector 

uj  Circular  frequency  of  the  applied  forces 

4>  Element  natural  frequency  mode  shape  (eigenvector) 
£  Element  structural  damping  coefficient 

Subscripts 
b  Baseline 

e  Element 

I  Imaginary  component  of  response 
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jth  element 

m  At  minimum  response 

n  At  resonance 

R  Real  component  of  response 

Superscripts 
T  Matrix  transpose 


